SCIENTIFIC 

REPpRTS 



»H5 




OPEN 



SUBJECT AREAS: 
GRAPHENE 
ATOMISTIC MODELS 

Received 
1 July 201 3 

Accepted 
1 August 201 3 

Published 
21 August 201 3 



Correspondence and 
requests for materials 
should be addressed to 
A.B. (angelo. 
bongiorno@chemistry. 

gatech.edu) 



Origin of the Chemical and Kinetic 
Stability of Graphene Oxide 



Si Zhou 12 & Angelo Bongiorno 1 



1 School of Chemistry and Biochemistry, Georgia Institute of Technology, 901 Atlantic Drive, Atlanta, Georgia, 30332-0400, USA, 
2 School of Physics, Georgia Institute of Technology, Atlanta 30332-0430, Georgia, USA. 

At moderate temperatures (< 70 °C), thermal reduction of graphene oxide is inefficient and after its 
synthesis the material enters in a metastable state. Here, first-principles and statistical calculations are used 
to investigate both the low-temperature processes leading to decomposition of graphene oxide and the role 
of ageing on the structure and stability of this material. Our study shows that the key factor underlying the 
stability of graphene oxide is the tendency of the oxygen functionalities to agglomerate and form highly 
oxidized domains surrounded by areas of pristine graphene. Within the agglomerates of functional groups, 
the primary decomposition reactions are hindered by both geometrical and energetic factors. The number of 
reacting sites is reduced by the occurrence of local order in the oxidized domains, and due to the close 
packing of the oxygen functionalities, the decomposition reactions become - on average - endothermic by 
more than 0.6 eV. 

Graphene oxide (GO) 1 ' 2 is a material with potential applications in ultra-thin electronics 3 ' 4 , optoelectro- 
nics 5 ' 6 , energy storage 7 " 10 , sensor devices 11 , and mechanical actuators 1213 . GO is generally obtained 
through the use of strong acid/base treatments 4 " 16 , resulting in the addition of oxygen in the forms of 
epoxide and hydroxyl species to the un-defected regions of the carbon layers, and carboxyl, carbonyl, and phenol 
groups to the edges of the carbon sheets 17 . As -synthesized GO is nonstoichiometric and hygroscopic, typically 
presenting 0:C ratios between 0.3 and 0.5, and an uncertain but variable amount of water molecules intercalated 
between the oxidized carbon layers 18 " 20 . Synthesized for the first time 150 years ago 21 , GO remains still today an 
elusive material system whose properties are difficult to control. 

Controlling the chemical structure of GO is crucial to, for instance, opening up a band gap for enabling 
applications in electronics 3 ' 4 , or to calibrate its catalytic 11 ' 22 and mechanical properties 1213,23 . Most of the efforts 
to control GO rely today on post-synthesis chemical and thermal treatments 24 " 31 , favoring the occurrence of 
reduction mechanisms, loss of oxygen content, and thus partial restoration of graphene-like properties. Thermal 
reduction of GO, in particular, initiates at temperatures larger than 70°C 30 , it leads to the desorption of H 2 0 and 
0 2 , but also CO and C0 2 30,32 , and consequently it is detrimental to the structural integrity of the graphene sheets 
and the physical properties of the reduced material. At temperatures lower than 70°C, reduction mechanisms 
become inefficient or die off, GO enters in a metastable state, and after its synthesis the material is driven towards 
a quasi- equilibrium state with a stable 0:C ratio and structure 33 ' 34 . While reduction mechanisms of GO and 
reduced GO have been investigated extensively in recent years 6,25 ' 35 " 38 , the stability and structure of GO aged or 
annealed at moderated temperature (< 70 °C) have received much less attention to date. Computations based on 
the density functional theory, in particular, have been used to generate and study model structures of GO 39 " 42 , to 
interpret experimental measurements 33 ' 34 ' 43 ' 44 , and to elucidate reduction mechanisms 38,42,45 and the thermodyn- 
amic stability of GO 46,47 . In spite of efforts, however, the origin of the chemical and kinetic stability of GO remains 
elusive to date. 

Here, we perform density functional theory (DFT) calculations, we develop an alternative and simplified 
energy scheme describing GO, and then we combine both types of calculations to elucidate the mechanisms 
underlying the chemical and kinetic stability of GO. In this work, we consider the case of graphene layers 
functionalized with epoxide and hydroxyl groups, and thus we investigate the stability of large- area GO 
sheets 33,34,48 . This study shows that the key factor underlying the stability of GO is the tendency of the oxygen 
functionalities to diffuse, agglomerate, and form highly oxidized domains surrounded by regions of pristine 
graphene. Within these agglomerates, spontaneous decomposition reactions producing 0 2 or H 2 0 molecules are 
hindered by both geometrical and energetic factors. In detail, the number of sites where pairs of oxygen functional 
groups can meet and react is reduced by the occurrence of local order in the highly oxidized domains, and due to 
the close packing of oxygen functional groups, at these sites the decomposition reactions become endothermic by 
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- on average - 0.6 eV. Overall, this study shows that both 0:C ratio 
and ageing processes can be used to enhance the chemical stability of 
GO, and therefore that aging as -synthesized GO at moderate tem- 
peratures may consist of a simple post- synthesis treatment to temper 
and control material properties, or produce new forms of GO. 

Results 

Figure 1 shows the possible elementary binary reactions between 
epoxide and hydroxyl species chemisorbed on a graphene layer lead- 
ing to desorption of 0 2 or H 2 0 molecules. The reaction energy dia- 
grams shown in Fig. 1 are derived from DFT calculations 49 . 
Transition states are calculated by using nudged elastic band 
(NEB) DFT calculations 49 . Technical details of the calculations (see 
Methods) and the full list of results are reported and discussed in the 
Electronic Supplementary Information (ESI). Throughout this work, 
free energy values of GO are referred to that ones of pristine graphene 
and a gas of 0 2 and H 2 0 molecules at room temperature and stand- 
ard pressure. In detail, the chemical potentials of 0 2 and H 2 0 are set 
equal to: 

Mo 2 /H 2 o = e dft [0 2 /H 2 O] + Afi 02 /H 2 o (!) 

where E DFT is the zero -temperature molecular energy computed 
from DFT, and A|i is an experimental correction accounting for 
the entropic and enthalpic contributions to the chemical potential 
of these gaseous molecular species 50,51 . The Gibbs free energy of 
formation, AG, of an arbitrary distribution of N e epoxides and N h 
hydroxyl groups on graphene is defined as: 

AG = E DFT [GO}-E DFT [G}-N e ^ -N h -03{^ +fi Hi0 } (2) 



where E DFT [GO] and E DFT [G] are the total energies computed by 
using our DFT scheme of defected and pristine graphene layers, 
respectively. Enthalpic and entropic contributions to AG associated 
to defect formation were estimated by using standard methods 51 in 
the case of a single epoxide species on graphene, a single hydroxyl 
species on graphene, and selected GO models consisting of disor- 
dered distributions of epoxide and hydroxyl species on graphene. 
These calculations show that, at room temperature and 1 atm, 
enthalpic and (vibrational) entropic corrections to AG in Eq. (2) 
amount to less than —0.1 eV per oxygen defect. For completeness, 
this latter energy value multiplied by the total number of oxygen 
functionalities is added to the right side of Eq. (2), thereby leading 
to an approximate but appropriate definition of AG. In Fig. 1 and 
throughout this work, values of AG are given per oxygen functional 
group. We note that the Gibbs free energies of formation of a single 
epoxide or hydroxyl group on graphene are found to be equal to 
about 1.1 eV and 1.6 eV, respectively (Fig. 1). These energy values 
indicate that sparse and homogeneous functionalizations of gra- 
phene with epoxide and hydroxyl groups are chemically unstable 
and prone to decompose via one of the binary reactions shown in 
Fig. 1 into 0 2 , H 2 0, and pristine graphene. 

The free-energy diagrams in Fig. 1 show that the most efficient - in 
energetic terms - low-temperature decomposition reactions of GO 
producing 0 2 and H 2 0 involve either two epoxides or two hydroxyls 
chemisorbed on the same side of the carbon basal plane. The reaction 
between the two epoxides is exothermic by more than 1 eV, it is 
controlled by an energy barrier of about 1 eV, and it involves an 
electronic singlet-to-triplet spin transition (ESI), while that one 
between the two hydroxyls is exothermic by about 0.5 eV and it 
involves energy barriers of about 0.5 eV (Fig. 1). The remaining 
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Figure 1 | Binary decomposition reactions between epoxide and hydroxyl groups on graphene. (a) Epoxide- epoxide reaction leading to formation of 0 2 . 
(b) Hydroxyl-hydroxyl reaction leading to formation of a H 2 0 molecule and an epoxide group on graphene. (c) Hydroxyl-epoxide reaction leading to 
formation of 0 2 and a hydrogen species chemisorbed on graphene. The reactions in (a), (b), and (c) involve oxygen species chemisorbed on the same side of 
the basal plane, (d) Two epoxides chemisorbed on opposite sides of a graphene layer forming a carbonyl-pair and then 0 2 . (e) Two hydroxyls on opposite 
sides of graphene transforming into a carbonyl-pair and hydrogen species, (f) An epoxide and a hydroxyl on opposite sides of the layer transforming into 
molecular oxygen and a hydrogen species. Energy values correspond to AG values per oxygen group computed via Eq. (2). Relevant energy barriers are 
computed by using NEB-DFT calculations. Bottom-right corner, legends of symbols used to represent reagents, by-products, and final products. 
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binary reactions between epoxides and hydroxyl groups lead to the 
formation of stable intermediate species that, with the exception of 
the carbonyl-pair group formed by the reaction of two epoxides 
chemisorbed on opposite sides of a graphene layer, lie higher in 
energy than the reagents by more than 1 eV. Carbonyl-pair groups 
are energetically more stable than two individual epoxide species on 
graphene (Fig. 1), and their transformation into 0 2 involves energy 
barriers as high as 4.9 eV, comparable to the energy costs required by 
oxygen functionalities to break through the graphene network and 
switch the side onto which they are chemisorbed. Thus, at moderate 
temperatures these reactions are unlikely to constitute viable path- 
ways leading to the decomposition of GO, which therefore is domi- 
nated by the reactions in Figs. 1(a) and 1(b). 

In addition to kinetic information, Figure 1 reports the energy 
gains yielded by associating two individual oxygen functionalities 
on graphene into a dimeric species formed by epoxide and/or 
hydroxyl groups sharing a C=C bond. The energy gains show, in 
particular, that the interaction between epoxide and hydroxyl species 
on graphene is, in all cases, attractive. The attraction is largest for 
hydroxyls, while for all possible combinations the association energy 
is larger when the two oxygen groups are chemisorbed on the oppos- 
ite rather than the same side of the carbon basal plane. The kinetics of 
agglomeration is governed by diffusion processes, and in case of 
isolated epoxide and hydroxyl groups on graphene, their diffusion 
rates are controlled by migration energies of 0.8 eV and 0.5 eV, 
respectively (ESI). A number of DFT studies have addressed the 
diffusion and interaction of epoxide and hydroxyl species on gra- 
phene, and the results thus far discussed are in agreement with these 
recent computational works 25 ' 37 ' 39 ' 40 ' 42 - 47 ' 52 - 56 . Here, we address in 
detail the effects played by agglomeration phenomena onto the 
chemical stability and structure of GO. 

Figure 2(a) shows the association energy of pairs of nearest neigh- 
bor epoxides, hydroxyls, and epoxide and hydroxyl groups chemi- 
sorbed on graphene. We consider, in particular, non-equivalent 
complexes formed by two oxygen functional groups, epoxide or 
hydroxyl species, separated by up to four C=C bonds and chemi- 
sorbed on either the same or opposite side of a basal plane. Our DFT 
calculations show that, in all cases, the two-species interaction is - 
on average - attractive and short-range, decaying to negligible 
values when the two oxygen groups are separated by more than 
4 C=C bonds (Fig. 2(a)). The association energies of pairs of oxygen 



functional groups separated by less than 4 C=C bonds are scattered 
over wide energy intervals (> 0.5 eV), indicating that the binding 
energy of such binary complexes results from a non-trivial competi- 
tion of strain and electronic effects. The two-species complexes and 
corresponding binding energies in Fig. 2(a) are here used to devise a 
simplistic and approximated energy scheme, alternative to DFT, and 
suited to describe arbitrary epoxide/hydroxyl functionalizations of 
graphene. 

We express the energy of a GO layer as sum of two-species energy 
contributions E IJK , where IJ refers to a pair of oxygen functional 
groups separated by less than 4 C=C bonds, and K refers to its 
corresponding binary association energy in Fig. 2(a). In detail, we 
write: 



AE< 2 > = 



IJ 



(3) 



where AE (2) is the total energy of the GO layer relative to the energy of 
individual oxygen functional groups on graphene. Based on Eq. (3), 
we can thus define the Gibbs free energy of formation of a GO layer 
hosting N e epoxides and N h hydroxyls as: 



AG® =AE® +N e n e + N h n h 



(4) 



where |i e and |i h are the Gibbs free energy of formation computed 
using DFT of a single epoxide and hydroxyl species on graphene, 
respectively. To test and calibrate this simple energy scheme, we use 
reference energies computed using DFT of an extended set of single- 
layer GO models, including model structures consisting of small 
agglomerates of a few epoxide or hydroxyl species on graphene to 
ordered and homogeneous structures of a GO layer presenting vary- 
ing 0:C ratios (ESI). The comparison between the energies obtained 
by using DFT and the simplified additive scheme shows that Eq. (3) 
reaches a good accuracy provided that, one, we exclude from the set 
of binary association energies reported in Fig. 2(a) that ones corres- 
ponding to epoxide-epoxide pairs exhibiting mild repulsive inter- 
action energies (< 0.3 eV), and two, we add a penalty energy of 
0.3 eV to each energy term E IJ)K associated to hydroxyl-hydroxyl 
complexes. With these two ad hoc corrections, our simplistic scheme 
reproduces energy values derived from DFT calculations of a variety 
of GO models with a mean absolute error of less than 0.1 eV 
(Fig. 2(b)). These comparisons also show that the aforementioned 
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Figure 2 | Pair-wise association energy of oxygen functional groups on graphene. (a) Energy computed using DFT (symbols) of binary epoxide -epoxide 
(red), hydroxyl-hydroxyl (blue), and epoxide-hydroxyl (cyan) complexes with functional groups separated by up to 4 C=C bonds. Energy values are 
obtained by referring the energy of the binary complex to that ones of pristine graphene and the isolated species on graphene. Insets show binary 
complexes having functional groups separated by one, two, and three C=C bonds (asterisks). Various non-equivalent binary complexes can be formed 
with oxygen functionalities separated by a selected number of C = C bonds (ESI). Solid colored lines are guide to the eye. (b) Energy of single-layer GO 
models presenting epoxy-only (red), hydroxyl-only (blue), and mixed epoxide-hydroxyl (cyan) functionalizations as computed from both DFT and the 
simplified energy scheme in Eq. (3). Values of AE (per functional group) computed using DFT are obtained by referring the energy of GO to that ones of 
pristine graphene and the isolated species on graphene. Filled symbols refer to crystalline structures with high 0:C ratios, while empty symbols show the 
energy of sparse agglomerates of functional groups on graphene with low 0:C ratios. Insets and ball and stick illustrations on the right side show selected 
GO structures considered in these comparisons; C, O, and H are shown in gray, red, and white colors (the full list of model structures is reported in the 
ESI). 
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ad hoc corrections are necessary to account for - approximately and 
on average - the decreasing and increasing strain energy of agglom- 
erates containing more than a pair of epoxide or hydroxyl groups, 
respectively (ESI). We use the aforementioned simplified energy 
scheme to predict the structure and energy of GO layers aged at mild 
temperatures (see Methods). 

Figures 3(a) and 3(b) report values of AG per functional group of 
epoxy- and hydroxyl -only functionalizations of graphene, respect- 
ively. Figure 3(a) reports Gibbs free energy values derived from Eq. 
(2), DFT calculations, and GO models presenting regular and homo- 
geneous distributions of oxygen functional groups (ESI), while 
Figure 3(c) reports energy values derived from Eq. (4), Monte 
Carlo simulated annealing simulations, and thus model structures 
of GO obtained by mimicking ageing at moderate temperatures (see 
Methods and ESI). Although greater than zero, the free energy values 
in Fig. 3(a) and 3(b) show that the chemical stability of GO in stand- 
ard ambient conditions increases with the 0:C ratio, nearly reaching 
equilibrium with the gas of 0 2 and H 2 0 molecules in ambient 
conditions in the case of fully oxidized and crystalline structures 
of GO. Furthermore, our Monte Carlo simulations show that, 
regardless the type of oxygen functionalization, ageing leads to the 
agglomeration of functional groups, enhanced chemical stability, 
and formation of a non-homogeneous phase of GO consisting of 
highly oxidized domains surrounded by areas of pristine graphene 
(Fig. 3). Similar results are obtained in the case of mixed epoxide - 
hydroxyl functionalizations of graphene (Fig. 3(c)). In this case, in 
addition to aggregation phenomena, ageing favors also segregation 
phenomena and formation within the highly oxidized domains of 
regions rich in epoxide and hydroxyl groups (Fig. 3(c)). DFT and 



lattice model Monte Carlo simulations therefore show that, due to 
the attractive interaction between oxygen functional groups on gra- 
phene, both 0:C ratio and ageing at mild temperature enhance 
chemical stability of GO and favor the appearance of partial struc- 
tural order. 



Discussion 

Our DFT and Monte Carlo calculations show that, in spite of the 
enhanced chemical stability, GO aged at moderate temperatures does 
not reach equilibrium with the gas phase and thus remains suscept- 
ible to decompose into gas molecules and pristine graphene (Fig. 3). 
Our calculations show, however, that aggregation phenomena and 
the formation of a non-homogeneous phase of oxidized and non- 
oxidized domains influence not only the chemical but also the kinetic 
stability of GO. In fact, the analysis of the models of aged GO gen- 
erated via Monte Carlo simulations shows that within the highly 
oxidized regions all the decomposition processes shown in Fig. 1 
become - on average - endo thermic, including the reactions between 
two nearest neighbor epoxides or hydroxyls chemisorbed on the 
same side of the basal plane; the energy diagrams of these two decom- 
position reactions are shown in Fig. 4. In these diagrams, the free 
energy of the initial state corresponds to AG of aged GO, that one of 
the final state is obtained by averaging the energy costs associated to 
decomposition reactions of eligible pairs of oxygen functionalities, 
while energy barriers are extrapolated from NEB -DFT calculations 
based on model structures of GO presenting a disordered distri- 
bution of epoxide and hydroxyl species (ESI). Our analyses - of a 
variety of model structures of aged GO generated by Monte Carlo 
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Figure 3 | Chemical stability and structure of large-area GO sheets, (a) AG per functional group of GO layers oxidized with only epoxide species. Black 
symbols connected by segments show values of AG computed using DFT and hand-made model structures of GO, encompassing small clusters and 
regular distributions of epoxide species of varying 0:C ratio. Red symbols connected by segments show values of AG of GO models generated from Monte 
Carlo simulated annealing simulations and thus mimicking aged structures of this material with varying 0:C ratios. Insets show a few selected model 
structures of aged GO functionalized with only epoxide groups. In these illustrations, graphene is not shown and red and yellow discs indicate epoxide 
species facing upward and downward the basal plane, respectively, (b) Same as (a) for hydroxyl-only functionalizations of graphene. In the insets, 
hydroxyl species facing upward and downward the carbon plane are shown in blue and cyan, respectively, (c) AG per functional group of GO layers 
oxidized with both epoxide and hydroxyl group. OH/O refers to the fraction of hydroxyl groups relative to the total amount of oxygen. The GO models are 
generated via Monte Carlo simulations (see Methods). Magenta and cyan symbols (connected by segments) refer to GO models with fixed 0:C ratios 
equal to 0.3 and 0.4, respectively. The illustrations on the right side show selected configurations of aged GO layers with an 0:C ratio of 0.4 and different 
relative fractions of hydroxyl and epoxide species. 
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Figure 4 | Binary decomposition reactions of aged GO. (a) Schematic energy diagram of the binary reaction between epoxide groups embedded in 
oxidized domains presenting high concentrations of epoxide groups. Values of AG (per functional group) of reagents and products correspond to 
averages extracted from the model structures of aged GO generated by Monte Carlo simulations and the approximated energy scheme in Eq. (3). The 
energy barrier of 1.6 eV shown in (a) refers to the energy value extracted from a NEB-DFT calculation using the model structure of GO on the right and 
the pair of reacting epoxide groups circled with the dashed magenta line, (b) Same as (a) in the case of two hydroxyls embedded in an oxidized region of 
aged GO layers rich in hydroxyl groups. Also in this case, values of AG for reagents and products are derived by using the energy scheme of Eq. (3) and 
model structures of aged GO. The energy barrier of 1.7 is extracted from a NEB-DFT calculation using the GO model shown on the right and the pair of 
reacting hydroxyl groups circled by the dashed cyan line. In the ball and stick illustration, the dashed box shows the planar dimensions of the supercell 
used in the DFT calculations. More details about this model structure and NEB-DFT calculations are reported in the ESI. 



simulated annealing simulations - show that within the highly 
oxidized regions, local order and spatial correlations reduce the 
fractions of epoxide and hydroxyl groups eligible to react to approxi- 
mately 0.3 and 0.2, respectively. The average heats of reaction of 
pairs of epoxides and hydroxyls forming 0 2 and H 2 0 are about 
0.6 eV and 0.8 eV, respectively. Further, our NEB-DFT calculations 
show that the energy barriers associated to the reactions shown in 
Fig. 4 attain, in the interior of oxidized domains, values larger than 
1.6 eV. These results concur in showing that the origin of the kinetic 
stability of graphene oxide stems from the attractive interaction 
between the oxygen functional groups. These interactions stabilize 
energetically the pairs of reacting species susceptible to form 0 2 or 
H 2 0, rendering the decomposition processes endothermic and thus 
unlikely. 

In summary, our first-principles calculations show that the most 
effective low-temperature decomposition reactions of GO are 
between pairs of epoxides or hydroxyls chemisorbed on the same 
side of graphene, that the interaction between oxygen functionalities 
on graphene is short-range and attractive, and that a simple and 
qualitatively correct energy scheme can be used to describe epoxy 
and/or hydroxyl functionalizations of graphene. First-principles 
and lattice model Monte Carlo simulations show that the key 
factor underlying the stability of GO is the tendency of the oxygen 
functionalities to diffuse, agglomerate, and form highly oxidized 
domains surrounded by regions of pristine graphene. Within 
these agglomerates, all the binary decomposition reactions pro- 
ducing 0 2 or H 2 0 molecules are hindered by both geometrical and 
energetic factors. In detail, the number of sites where pairs of oxygen 
functional groups can meet and react is reduced by the occurrence of 
local order in the highly oxidized domains, and due to the close 
packing of oxygen functional groups, at these sites any decomposi- 
tion reaction becomes endothermic by - on average - more than 
0.6 eV. Overall, this study shows that both 0:C ratio and ageing 
processes can be used to enhance both chemical and kinetic stability 
of GO. Aging as -synthesized GO at moderate temperatures may 
therefore consist of a simple post- synthesis treatment to both temper 



and control material properties, as well as to lead to new forms of 
GO. 

Methods 

Density functional theory calculations are performed by using the QUANTUM- 
Espresso software 49 . We use an energy cutoff of 70 Ry, norm-conserving pseudopo- 
tentials 57 for all atomic species, and the parametrization of the exchange and 
correlation energy functional proposed by Perdew, Burke, and Ernzerhof 58 . Test 
calculations show that cutoff energy, k-point mesh in the Brillouin zone, and size and 
geometry of supercells introduce errors of about ±0.1 eV on our energy values. Bond 
energies of 0 2 , CO, and H 2 0 given by this first principles scheme are 5.3 eV, 10.9 eV, 
and 4.8 eV, respectively, within 3% of experimental data. First principles calculations 
of crystalline structures of GO are based on the use of small unit cells and dense k- 
point meshes, while the rest of the calculations are carried out by using supercells with 
planar dimension equal or larger than 5X6 graphene unit cells and by sampling the 
Brilloing zone in the T-point. The vacuum region separating the periodic graphene 
layers is - in all cases - equal to 12 A. 

Binding energy values and geometries of the binary complexes used to calculate 
via Eq. (3) the energy of an arbitrary distribution of epoxide and hydroxyl species 
chemisorbed on graphene are reported in the ESI. This simplified energy scheme 
is used in conjunction with lattice model representations of functionalized gra- 
phene and kinetic Monte Carlo simulations. In this numerical framework, epox- 
ide- and hydroxyl-like defects occupy C=C bonds and C sites on one of the two 
opposite sides of the graphene plane, and Eq. (3) is used to compute the energy of 
a distribution of oxygen defects. The functional groups can move freely - with no 
energy barriers - between nearest neighbor defect sites. A Monte Carlo step 
consists in displacing a functional group selected at random to an empty nearest 
neighbor defect site, this one also selected at random. The energies computed by 
using Eq. (3) of the initial and modified defect distributions are used to calculate a 
Boltzmann factor, and thus accept or reject the new configuration. This lattice 
model Monte Carlo scheme is used in conjunction with simulated annealing 
simulations to mimic the ageing of GO and predict the structure resulting from 
the agglomeration of homogeneous distributions of functional groups (ESI). These 
Monte Carlo simulations are carried out by using large models of GO, including 
hundreds of oxygen functionalities. Number and nature of the oxygen functional 
groups are held fixed in the simulations. 

It is to be noted that in our Monte Carlo simulations, temperature is a 
numerical parameter used to control (accelerate and/or arrest) the occurrence of 
structural transformations in our GO models. Temperature in our Monte Carlo 
simulations is thus not intended to mimic any experimental temperature. Due to 
the nature of the underlying energy scheme - which accounts only for the 
interaction of hydroxyl and epoxide species on a graphene layer, raising the 
temperature within our conventional Monte Carlo scheme does not introduce any 
new physics or chemistry besides that one described by the discrete energy 
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scheme, and the simulated annealing simulations correspond simply to a 
numerical strategy used to speed up structural transformations and determine - 
by a slow quenching - a low free-energy arrested state of GO. Therefore, our 
Monte Carlo simulations have no direct link with the experiment; they are simply 
used to speed up the evolution (ageing) of GO in a regime where structural 
transformations are driven only by the diffusion of the epoxides and hydroxyls, i.e. 
moderate (experimental) temperatures. 



1. Dikin, D. A. et al. Preparation and characterization of graphene oxide paper. 
Nature (London) 448, 457-460 (2007). 

2. Dreyer, D. R., Park, S., Bielawski, C. W. & Ruoff, R. S. The chemistry of graphene 
oxide. Chem. Soc. Rev. 39, 228-240 (2010). 

3. Wei, Z. et al. Nanoscale tunable reduction of graphene oxide for graphene 
electronics. Science 328, 1373-1376 (2010). 

4. Wu, X. et al. Epitaxial-graphene/graphene-oxide junction: an essential step 
towards epitaxial graphene electronics. Phys. Rev. Lett. 101, 026801 (2008). 

5. Loh, K. P., Bao, Q., Eda, G. & Chhowalla, M. Graphene oxide as a chemically 
tunable platform for optical applications. Nature Chem. 2, 1015 (2010). 

6. Kumar, P. V., Bernardi, M. & Grossman, J. C. The Impact of Functionalization on 
the Stability, Work Function, and Photoluminescence of Reduced Graphene 
Oxide. ACS Nano 7, 1638-1645 (2013). 

7. Chen, Y., Zhang, X., Zhang, D., Yu, P. & Ma, Y. High performance supercapacitors 
based on reduced graphene oxide in aqueous and ionic liquid electrolytes. Carbon 
49,573-580 (2011). 

8. Gao, W. et al. Direct laser writing of micro- supercapacitors on hydrated graphite 
oxide films. Nature Nanotechnology 6, 496 (2011). 

9. Stoller, M. D., Park, S. J., Zhu, Y. W. & Ruoff, R. S. Graphene-Based 
Ultracapacitors. Nano Lett. 8, 3498-3502 (2008). 

10. Zhu, Y. et al. Carbon-Based Supercapacitors Produced by Activation of Graphene. 
Science 332, 1537-1541 (2011). 

11. Robinson, J. T., Perkins, F. K., Snow, E. S., Wei, Z. Q. & Sheehan, P. E. Reduced 
Graphene Oxide Molecular Sensors. Nano Lett. 8, 3137 (2008). 

12. Park, S. et al. Graphene Oxide Papers Modified by Divalent Ions — Enhancing 
Mechanical Properties via Chemical Cross-Linking. ACS Nano 2, 572-578 
(2008). 

13. Suk, J. W., Piner, R. D., An, J. H. & Ruoff, R. S. Mechanical Properties of Mono 
layer Graphene Oxide. ACS Nano 4, 6557-6564 (2010). 

14. Hummers, W. S. & Offeman, R. E. Preparation of Graphitic Oxide. J Am Chem Soc 
80, 1339-1339 (1958). 

15. Dimiev, A., Kosynkin, D. V., Alemany, L. B., Chaguine, P. & Tour, J. M. Pristine 
Graphite Oxide. /. Am. Chem. Soc. 134, 2815-2822 (2012). 

16. Marcano, D. C. et al. Improved Synthesis of Graphene Oxide. ACS Nano 4, 
4806-4814 (2010). 

17. Lerf, A., He, H., Forster, M. & Klinowski, J. Structure of Graphite Oxide Revisited. 
/. Phys. Chem. B 102, 4477-4482 (1998). 

18. Cerveny, S., Barroso-Bujans, F., Alegria, A. & Colmenero, J. Dynamics of Water 
Intercalated in Graphite Oxide. /. Phys. Chem. C 114, 2604-2612 (2010). 

19. Buchsteiner, A., Lerf, A. & Pieper, J. Water dynamics in graphite oxide 
investigated with neutron scattering. / Phys Chem B 110, 22328-22338 (2006). 

20. Lerf, A. et al. Hydration behavior and dynamics of water molecules in graphite 
oxide. / Phys Chem Solids 67, 1 1 06- 1 1 1 0 (2006) . 

21. Brodie, B. C. On the Atomic Weight of Graphite. Phil. Trans. Royal Soc. London A 
149, 249-259 (1859). 

22. Dreyer, D. R., Jia, H. & Bielawski, C. W. Graphene oxide: a convenient 
carbocatalyst for facilitating oxidation and hydration reactions. Angewl Cheml 
Intl Ed. 49, 6813-6816 (2010). 

23. Robinson, J. T. et al. Wafer- scale Reduced Graphene Oxide Films for 
Nanomechanical Devices. Nano Lett. 8, 3441 (2008). 

24. Acik, M. & Chabal, Y. J. Reducing graphene oxide for band gap engineering. Appl 
SurfSci, Submitted (2011). 

25. Bagri, A. et al. Structural evolution during the reduction of chemically derived 
graphene oxide. Nature Chem. 2, 581-587 (2010). 

26. Eda, G., Fanchini, G. & Chhowalla, M. Large-area ultrathin films of reduced 
graphene oxide as a transparent and flexible electronic material. Nature 
Nanotechnology 3, 270-274 (2008). 

27. Ekiz, O. 6., Urel, M., Giiner, H., Mizrak, A. K. & Dana, A. Reversible electrical 
reduction and oxidation of graphene oxide. ACS Nano 5, 2475-2482 (2011). 

28. Gomez-Navarro, C. et al. Electronic Transport Properties of Individual 
Chemically Reduced Graphene Oxide Sheets. Nano Lett. 7, 3499-3503 (2007). 

29. Gong, C, Acik, M., Abolfath, R. M., Chabal, Y. & Cho, K. Graphitization of 
Graphene Oxide with Ethanol during Thermal Reduction. The Journal of Physical 
Chemistry C 116, 9969-9979 (2012). 

30. Jung, I. et al. Reduction kinetics of graphene oxide determined by electrical 
transport measurements and temperature programmed desorption. /. Phys. 
Chem. C 113, 18480-18486 (2009). 

31. Hossain, M. Z. et al. Chemically homogeneous and thermally reversible oxidation 
of epitaxial graphene. Nat. Chem. 4, 305-309 (2012). 

32. Yang, D. et al. Chemical analysis of graphene oxide films after heat and chemical 
treatments by X-ray photoelectron and Micro-Raman spectroscopy. Carbon 47, 
145-152 (2009). 



33. Kim, S. et al. Room-temperature metastability of multilayer graphene oxide films. 
Nature Materials 11, 544-549 (2012). 

34. Zhou, S., Kim, S. & Bongiorno, A. Chemical Structure of Oxidized Multilayer 
Epitaxial Graphene: A Density Functional Theory Study. The Journal of Physical 
Chemistry C 117, 6267-6274 (2013). 

35. Acik, M. et al. Unusual infrared-absorption mechanism in thermally reduced 
graphene oxide. Nature Materials 9, 840-845 (2010). 

36. Gomez-Navarro, C. et al. Atomic structure of reduced graphene oxide. Nano Lett. 
10, 1144-1148 (2010). 

37. Ghaderi, N. & Peressi, M. First- Principle Study of Hydroxyl Functional Groups on 
Pristine, Defected Graphene, and Graphene Epoxide. J Phys Chem C 114, 
21625-21630 (2010). 

38. Larciprete, R., Fabris, S., Sun, T. P., Lacovig, A. B. & Lizzit, S. Dual path 
mechanism in the thermal reduction of Graphene oxide. /. Am. Chem. Soc. 133, 
17315-17321 (2011). 

39. Yan, J.-A. & Chou, M. Y. Oxidation functional groups on graphene: structural and 
electronic properties. Phys. Rev. B 82, 125403 (2010). 

40. Yan, J.-A., Xian, L. & Chou, M. Y. Structural and Electronic Properties of Oxidized 
Graphene. Phys. Rev. Lett 103, 086802 (2009). 

41. Johari, P. & Shenoy, V. B. Modulating Optical Properties of Graphene Oxide: Role 
of Prominent Functional Groups. ACS Nano 5, 7640-7647 (2011). 

42. Paci, J. T., Belytschko, T. & Schatz, G. C. Computational Studies of the Structure, 
Behavior upon Heating, and Mechanical Properties of Graphite Oxide. /. Phys. 
Chem. Clll, 18099-18111 (2007). 

43. Kudin, K. N. et al. Raman spectra of graphite oxide and functionalized graphene 
sheets. Nano Letters 8, 36-41 (2008). 

44. Lu, N., Huang, Y., Li, H. B., Li, Z. Y. & Yang, J. L. First principles nuclear magnetic 
resonance signatures of graphene oxide. Journal of Chemical Physics 133, 034502 
(2010). 

45. Sun, T., Fabris, S. & Baroni, S. Surface Precursors and Reaction Mechanisms for 
the Thermal Reduction of Graphene Basal Surfaces Oxidized by Atomic Oxygen. 
Jl Physl Cheml C 115, 4730-4737 (2011). 

46. Wang, L. et al. Stability of graphene oxide phases from first-principles 
calculations. Phys. Rev. B 82, 161406 (2010). 

47. Liu, L. Z. et al. Amorphous structural models for graphene oxides. Carbon 50, 
1690-1698 (2012). 

48. Zhao, J. P., Pei, S. F., Ren, W. C, Gao, L. B. & Cheng, H. M. Efficient Preparation of 
Large-Area Graphene Oxide Sheets for Transparent Conductive Films. ACS Nano 
4, 5245-5252 (2010). 

49. Giannozzi, P. etal. QUANTUM ESPRESSO: a modular and open-source software 
project for quantum simulations of materials. /. Phys.: Condens. Matter 21, 395502 
(2009). 

50. Codata Recommended Key Values for Thermodynamics, 1977 - Report of 
Codata-Task-Group on Key Values for Thermodynamics. 1977. / Chem 
Thermodyn 10, 903-906 (1978). 

51. Reuter, K. & Scheffler, M. Composition, structure, and stability of Ru02(l 10) as a 
function of oxygen pressure. Phys Rev B 65, 035406 (2002). 

52. Boukhvalov, D. W. & Katsnelson, M. I. Modeling of Graphite Oxide. /. Am. Chem. 
Soc. 130, 10697-10701 (2008). 

53. Sljivancanin, Z., Milosevic, A. S., Popovic, Z. S. & Vukajlovic, F. R. Binding of 
atomic oxygen on graphene from small epoxy clusters to a fully oxidized surface. 
Carbon 54, 482-488 (2013). 

54. Nguyen, M. T., Erni, R. & Passerone, D. Two-dimensional nucleation and 
growth mechanism explaining graphene oxide structures. Phys Rev B 86, 
115406 (2012). 

55. Huang, B., Xiang, H. J., Xu, Q. & Wei, S. H. Overcoming the Phase Inhomogeneity 
in Chemically Functionalized Graphene: The Case of Graphene Oxides. Phys Rev 
Lett 110, 085501 (2013). 

56. Topsakal, M. & Ciraci, S. Domain formation on oxidized graphene. Phys Rev B 86, 
205402 (2012). 

57. Troullier, N. & Martins, J. L. Efficient pseudopotentials for plane-wave 
calculations. Phys. Rev. B 43, 1993-2006 (1991). 

58. Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized Gradient Approximation 
Made Simple. Phys. Rev. Lett. 77, 3865-3868 (1996). 



Acknowledgements 

The authors acknowledge the support by the National Science Foundation (NSF) through 
the grants CMMI-1 100290, DMR-0820382, and CHE-0946869, and the partial support by 
the Samsung Advanced Institute of Technology (SAIT). The authors also thank Suenne 
Kim, Elisa Riedo, Claire Berger, Walt de Heer, and Yves J. Chabal for the frequent and useful 
discussions. 

Author contributions 

S.Z. carried out calculations and analyses. A.B. supervised the research work and wrote the 
manuscript. Both authors reviewed the manuscript. 



SCIENTIFIC REPORTS | 3 : 2484 | DOI: 1 0.1 038/srep02484 



6 



Additional information 

Supplementary information accompanies this paper at http://www.nature.com/ 
scientificreports 

Competing financial interests: The authors declare no competing financial interests. 



How to cite this article: Zhou, S. & Bongiorno, A. Origin of the Chemical and Kinetic 
Stability of Graphene Oxide. Set Rep. 3, 2484; DOI:10.1038/srep02484 (2013). 

l/££v ^ I This work is licensed under a Creative Commons Attribution 3.0 Unported license. 
K-^^KZ^H To view a copy of this license, visit http://creativecommons.Org/licenses/by/3.0 



SCIENTIFICREPORTS | 3 : 2484 | DOI: 1 0.1 038/srep02484 



7 



